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A narrow  band  wave  spectrum  was  applied  to  theoretical 
relationships  previously  developed  for  set  down  and  set  up 
in  an  attemp t to f i n d second  order  non-steady  solutions 
for  these  concepts.  The  initial  effort  was  to  apply  this 
spectrum  to  the  radiation  stress  tensor  using  linear  wave 
theory.  Another  development  was  attempted  by  i ncoroorat  i ng 
the  spectrum  and  the  solution  of  the  long  wave  equation  in- 
to the  Bernoulli  and  vertical  momentum  equations.  Results 
obtained  indicate  that  the  solution  for  mean  water  level 
outside  the  surf  zone  is  composed  of  a steady  component  and 
a periodic  unsteady  component;  the  periodic  component  being 
of  fhe  form  of  a long  wave  with  a frequency  lower  tnan  the 
components  of  the  wave  spectrum.  The  solution  for  set  up 
is  then  composed  of  the  same  type  components.  The  exact 
relationships  depend  on  the  patching  process  that  is  made 
for  the  solutions  through  the  breaker  line. 
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I NTRODUCTI ON 


Mos+  of  the  observable  phenomena  along  a coastline  are 
the  direct  result  of  the  action  of  the  incoming  waves, 
waves  which  begin  in  most  instances  as  a disorganized  con- 
fused state  of  the  ocean  surface,  produced  by  a storm  far 
at  sea.  In  their  transit  across  the  vast  expanse  of  the 
ocean,  they  begin  to  sort  themselves  out  and  form  a somewhat 
regular  oscillation  of  the  ocean  surface,  the  lower  freauen- 
cy  oscillations  traveling  faster  and  thereby  leading  the 
train.  This  train  eventual  ly  ends  by  encountering  a beacn, 
where  its  energy  is  expended  in  the  form  of  breaking  waves. 

Since  it  is  this  aspect  of  the  wave's  life  cycle  that 
influences  man  the  most,  consideraole  effort  has  been  ex- 
pended investigating  this  area.  Among  other  effects,  it  nas 
been  observed  that  waves,  in  tne  process  of  shoal  ing  and 
eventual  breaking,  produce  a variation  in  The  mean  sea  level 
This  variation  in  sea  level  has  been  considered  as  the  pri- 
mary cause  for  sucn  nearsnore  currents  as  rip  currents.  The 
sea  surface  variation  consists  of:  (a)  a graaual  depression 
of  the  mean  sea  level  beginning  offshore  and  reaching  a maxi 
mum  a+  the  breaker  line  and  (b)  inside  the  surf  zone  a slope 
of  mean  sea  level  which  increases  and  extends  shoreward  to 
a point  on  the  beach  higher  than  the  still  water  line.  The 
depression  is  termed  set  down  and  the  slope  is  called  set  up 


It  is  reasonable  to  expect  that  variations  in  the  amount  of 
set  down  or  set  up  along  a beach  can  provide  the  head  to 
produce  a current. 

Previous  investigations  have  considered  steady  s'^ate 
solutions.  Experimental  results  obtained  by  both  3owen 
D967^  and  Van  Dorn  \i\916^  agree  quite  favorably  with  the 
steady  state  solutions  produced  by  Longuet-H i gg  i n s and 
Stewart  [I962j  using  linear  wave  theory. 

This  investigation  considers  an  application  of  a simple 
wave  spectrum  to  the  existing  theories  in  an  attempt  to  ob- 
tain a non-steady  solution  for  the  set  down  and  set  up 
phenomena.  An  opening  chapter  on  baokground  is  provided  to 
ensure  the  necessary  understanding  of  the  existing  theories 
and  their  development.  Included  is  a section  on  the  develop- 
ment of  the  "radiation  stress  tensor",  a concept  which  was 
proved  useful  by  Longuet-H i gg i ns  and  Stewart  Ll962j  in 
treating  the  shoaling  process  of  waves. 

Chapter  III  deals  with  the  application  of  a simple 
spectrum  to  the  derivation  of  set  down  and  set  ud.  First 
order  linear  theory  and  first  order  tneorv  Including  a sloa- 
ing  bottom  are  used  to  describe  the  spectral  wave  comconents. 
The  final  chap+ers  conclude  with  a comparison  or  the  steady 
state  solutions  produced  by  the  earl ier  work  and  the  unsteaoy 
results  obtained  here.  The  numerical  results  are  comoared 
with  tne  results  for  the  steady  case  given  oy  Bowen  iIl96T]. 


1 I . BAlKGROUND 


Changes  in  mean  sea  level  near  a.  shoreline  have  been 
studied  both  theoretically  and  experimentally.  A theoreti- 
cal framework  was  formulated  by  Long uet-H i gg ! n s and  Stewart 
ri962,  1963,  19643  which  dealt  with  the  excess  momentum 

flux  due  to  the  presence  of  unsteady  wave  motion  and  which 
they  termed  "radiation  stress".  Longuet-H i gg i n s and 
Stewart  were  able  to  define  many  of  the  shoaling  effects  of 
a train  of  waves  including  wave  set  down  and  set  uo  utiliz- 
ing this  radiation  stress  concept.  This  chapter  reviews 
the  development  of  the  radiation  stress  tensor  and  its  re- 
lationship to  the  conceots  of  set  down  and  set  up.  It  also 
includes  a direct  approach  used  by  Longue+-Higgins  315673  to 
derive  an  expression  for  set  down  utilizing  the  vertical 
momentum  equation  and  the  Bernoulli  equation. 


A.  DEVELOPMENT  OF  THE  RADIATION  STRESS  TENSOR 

Since  the  approach  of  Lo n g u e t-K  i gg  i n s and  Stewart  is 
rather  lengthy  and  tends  to  obscure  the  conceots  involved, 
the  later  and  more  systematic  approach  by  Phi  I I ios  315663  is 
used.  The  development  of  the  radiation  stress  tensor  and 
*ne  resulting  phenomena  of  set  down  and  set  up  is  kep*  as 
general  as  oossible.  The  coordinate  system  is  given  by 
Figure  (1)  where  the  x-axis  is  p e rpen  d i c u 1 a r to  '^he  snore, 
the  y-axis  is  parallel  to  the  snore,  and  tne  z-axis  is 
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vertically  upward  from  the  still  water  level.  The  govern- 
ing equations  are  the  continuity  equation  and  the  hori- 
zontal momentum  equations. 

The  continuity  equation  is  given  by 


3p 

3t 


+ 


apu . 

3x  . 

I 


3pw 

IT-  = 


0 , 


i = I,  2 


( I ) 


where  I,  2 refer  to  the  x,  y components.  The  horizontal 
velocity  u.  is  composed  of  a mean  flow  component,  Uj,  and 
a fluctuating  component  representing  the  deviation  from 
mean  flow,  u.',  such  that  u . = U . + u.'.  Since  there  is 
no  mean  flow  in  the  vertical,  w = w'.  Multiplying  the  con- 
tinuity equation  by  u.  and  adding  this  result  to  the  hori- 
zontal momentum  equation 


3 u . 

1 

3u  . 

3 u . 

1 

3t 

+ PU  . -r + 

J 

3 0 

3x  . 
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+ R. 
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(2) 


p rod  uces 


3u  . 

I 

3t 


3 p u . u . 3 p u . w 

' J-  + 


9x  . 
J 


3 z 


- t R. 

3x . I 


( 3 ) 


Integrating  over  depth  from  -h  to  q , using  Leibnitz's  rule 
and  applying  the  kinematic  free  surface  and  bottom  boundary 
condition  yields 
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By  time  averaging  this  equation  term  by  term  and  by  making 
the  following  definitions 


M.  = / pU.dz  = pU.(ri  + h)  , 


-h 


pu  . ’ dz  , 


M.  = M.  + M.  , 

I I I 


M . M . 

U . = ^ = U . + ! 

' p(r|  + h)  ‘ p(fi'+h) 


the  following  expression  is  obtained: 


3M.  » ^ ^ r M.M. 

■5^  + 4—  {U.M.  + / LPu.'u.'  + p5..Ildz  - '■ — ^ 

J J p(n  + h) 

-h 

- J pg  (n  + h ) 3 . } 


T,  . R,  . 


5..=0for  i 

= I for  i = j 


The  first  term  on  the  left  hand  side  is  the  local  change  in 
horizontal  momentum  flux.  In  the  braces,  the  first  term 


represents  the  momentum  flux  produced  by  the  steady  state 
flow.  The  last  three  terms  in  the  bra'us  contain  all  the 
unsteady  contribution  to  the  momenum  flux  with  -f-he  Hydro- 


static effect  subtracted  out;  this  is  the  momentum  flux  due 
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to  the  unsteady  motion  or  the  excess  momentum  flux  referred 
to  as  the  radiation  stress  tensor: 


J Cpu. 'u^. '+p6,^.]d. 


M.M  . 


p ( n+h ) 


J pg  ( n + h ^ j j 


Equation  (5)  is  then  simplified  to 


It-M.  +1-  {U.M.  +S..}  = T.  +R. 

at  I ax^.  I I I j I I 


T.  represents  the  horizontal  force  produced  by  the  slope  of 
the  free  surface  and  is  given  by 


Ti  = - pg(n+h)  . 


Outside  the  surf  zone  it  is  assumed  n <<  -h  and  (8)  reduces 


"i  = - • 


R.  is  the  averaged  and  integrated  frictional  stress  term. 
These  equations  are  general  ana  aoply  to  all  kinds  of  steady 
and  unsteady  motion.  The  only  simplifying  assumption  is 
that  mean  flow  is  uniform  over  depth.  The  advantage  of  us- 
ing the  "radiation  stress"  technique  to  solve  physical  prob- 
lems is  that  the  second  order  effects  are  obtained  using 
first  order  theory. 


B.  WAVE  SET  DOWN 

To  present  the  concept  of  set  down,  consider  only  the  x- 
component  of  horizontal  momentum  flux  equation  (7)  with 
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waves  that  propagate  shoreward  from  deep  water  with  their 
crests  making  an  arbitrary  angle  ct  with  the  shoreline.  For 
simplicity,  the  bottom  is  assumed  to  be  composed  of  parallel 
contours  so  that  gradients  in  the  y-direction  are  zero.  The 
bottom  slope  is  allowed  to  vary  only  gradually  so  that 
energy  reflection  from  the  shore  may  be  neglected  and  the 
shoaling  effects  caused  by  changes  in  the  bottom  can  be  con- 
sidered in  a step-like  fashion.  By  neglecting  the  frictional 
effects  and  assuming  that  any  current  gradient  in  the  x- 
direction  is  small,  (7)  can  be  written  for  outside  the  surf 
zone  as 


9S 


XX 


3x 


(10) 


Since  (10)  is  an  equation  in  two  unknowns  (M^  and  n ) > 
a second  equation  is  required  for  a solution.  To  provide 
it,  the  continuity  equation,  (I),  is  vertically  integrated 
over  depth  by  the  use  of  Leibnitz's  rule.  This  produces  the 
conservation  of  mass  flux  equation  as  given  by  Phillips 
Cl966]]  where  gradients  in  the  y-direction  are  zero. 
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9n 

TF 


+ 


aM 

X 

ax 


0 . 


(II) 


utilizing  linear  wave  theory,  the  radiation  stress  tensor 
is  proportional  to  the  sauare  of  the  local  wave  amplitude,  a^ 
or  to  the  total  energy.  Therefore,  following  rhe  method  of 

the 


Long  ue  t-H  i gg  i n s and  S + ewart  ill962j,  the  applied  force  of 


system  travels  with  the  speed  of  the  energy,  i.e.. 


at  the 


group  velocity,  c^.  Hence  the  transformation  9/9t  + 

Cg  9/9x  = 0 can  be  applied.  Equations  (10)  and  (II)  become 


X , u on 

^ -97 


9S 


XX 


dX 


X or) 

9x  ^g  Tx 


for  which  the  solution  is 


= 0 , 


(12) 


(13) 


In 

9x 


or  on  integration 


9S 


^ ( gh  - ) 


XX 


9x 


(14) 


" = P 


, , 2 XX 

( g h - c ) 

9 


(15) 


Since  Cg  ->■  gh  in  shallow  water,  equation  (15)  is  a non- 
steady state  solution  which  implies  that  the  mean  water 
level  increases  negatively  withou-^  bound  as  the  wave  moves 
into  shallow  water.  The  explanation  offeree  by  Longuer- 
Higgins  and  Stewart  [I952j  for  this  apparent  resonant  condi- 
tion is  that  its  effect  takes  time  to  build  and  the  energy 
involved  Is  dissipated  prior  to  it  reaching  any  significance 
by  the  breaking  of  the  wave. 

By  imposing  steady  state  conditions  and  descriping  One 
unsteady  motion  using  linear  wave  theory,  Longuet-H i gg i ns 
and  Stewart  Cl962j  subsequently  developed  the  following  si~u- 
at  i on  for  n : 


I 


n = - T 


a k 


2 s i n h 2 kh 
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(16) 


This  is  a second  order  equation  in  terms  of  the  local  depth, 
amplitude,  and  wave  number.  It  is  apparent  that  as  the  wave 
train  approaches  the  point  of  breaking,  ii  decreases  and  there 
is  a set  down  of  the  mean  water  level.  The  set  down  in- 
creases to  the  point  where  the  wave  breaks  and  other  assump- 
tions regarding  n must  be  applied. 

Longuet-H  i gg  i ns  [II9673  also  derived  this  solution  in 
another  manner  without  referring  to  the  radiation  stress 
term.  By  integrating  the  vertical  momentum  equation  over 
depth  and  time  averaging,  the  total  average  vertical  momentum 


is  obtained: 

( p - pw^ ) 2_q  - pgn  = 0 . (17) 

A second  equation  used  is  Bernoulli's  integral 

p + i p(u  ^ + u ^ + w^)  + pgz  + p -U  = 0 , (13) 

2x  y 3t  ' 

with  the  restriction  That  the  flow  be  i r rota t i ona  I . By  set- 
ting z = 0 and  time  averaging,  (13)  becomes 

p - -^ip(u^  + n^  + w^)  - + C = 0,  (19) 

'^z  = 0 2 X y z = 0 ' 

where  C is  at  most  a constant.  From  (17)  and  (19),  ? ,, 

2 = J 

can  be  eliminated  giving 

p=_J-i(u2  + u2_  + C = 0 . (20) 

“ d X y 2 - u 
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From  this,  the  difference  in  mean  sea  level  is  obtained  at 
two  different  points  (Xj,  0)  and  (x^,  y ^ thereby 

eliminating  the  constant  C, 


AS  = - jL  ^ -^2  . , (21, 


The  velocities  are  expressed  applying  linear  wave  theory  us- 
ing the  following  relationships: 


^ K ( z-h  ) cos  ( kx  ' -at-!-e  ) , 

X s 1 n h kh  ' 


- ^ ~rr  cosh  k ( z-h  ) cos  ( kx  ' -at  + e ) ,(22) 

s I n h kh 


s i n h kh 


sinh  k ( z- h ) s i n ( kx ' -at+e ) , 


where  x'  = x cosa  + y s i na  and  e denotes  an  arbitrary  con- 
stant phase  angle.  Substituting  equations  (22)  Into  the 
right  hand  side  of  (21)  yields 


I .-  2 - 2 -2. 

■= — u + u - w ; = 
2g  X y 


2„2 
a a 


4g  s i n h “ k h 


and  the  difference  in  mean  sea  level  is  then  obtained  as 


Ar  - a^k  ^ I 

" ■ 2s  i nh  2kh-'2  * 
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By  assuming  the  point  (X2,  'i  2*  ‘ infinitely  deep 

water  so  that  r\  2 = 0,  the  value  for  n at  the  point  .x^,  y ^ , 


0)  becomes 


2sinh2kh' 


This  solution  is  very  simple  and  straightforward,  with  the 
same  assumptions  as  before  applying  to  the  bottom  slope 
and  local  depth  h.  This  solution  is  plotted  In  Figure  (2) 
with  the  strong  correlation  to  laboratory  results  given  by 
Bowen,  inn man,  and  Simmons  [19633. 

As  waves  progress  shoreward  they  travel  from  deep  water 
to  shallow  water.  The  transition  in  the  linear  equation 
can  be  made  by  utilizing  the  conservation  of  energy  flux 
CE  • Cg  = constant).  However,  when  a wave  approaches  +he 
point  of  breaking, the  wave  steepens  beyond  that  allowed  by 
linear  theory  and  another  approximation  must  be  made.  It 
is  assumed  that  at  breaking 


where  y is  a proportionality  constant  which  is  generally 
assumed  to  be  equal  to  0.73  borrowed  from  solitary  wave 
theory.  Then  the  mean  sea  level  decression  at  "I'he  breaking 
point,  is  given  by 


^ - 76  "b  • 
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SET  UP 


Set  up  is  a phenomenon  that  occurs  in  the  surf  zona 
shoreward  of  the  breaker  line.  Because  energy  decreases 
shoreward  in  this  zone,  a different  formulation  for  the 
wave  amplitude  must  be  used  than  the  local  value  used  be- 
fore. It  is  assumed  that  the  breakers  being  considered 
are  of  the  spilling  type  which  retain  their  harmonic  charac- 
teristics and  gradually  decrease  in  amplitude  as  they 
progress  shoreward  such  that  wave  height  is  defined  by 

H = (28) 


k 


By  assuming  steady  state  conditions,  (29)  can  be  written  as 


i_  5 E = 

3x  2 


- pg(n  + h)  . 


Since  energy  decreases  shoreward,  it  is  evident  that  n in- 
creases and  produces  a set  up  effect.  Solving  for  n 


n = - kh  + C , 


where 


3 2 
8 ^ 


By  applying  the  value  of  ri|^  at  the  breaker  line,  the  con- 
stant C can  be  el  iminated  and 


ri  = k(h,-  h)  + nK  . 

0 D 


Where  a beach  has  a constant  slope  given  by  h = mx,  the  set 
up  also  has  a constant  slope  proportional  to  the  beach  slope. 
It  should  be  noted  that  the  simple  solution  obtained  by 
application  of  Bernoulli's  integral  oulside  the  surf  zone 
is  not  applicable  since  the  motion  is  no  longer  irrotational 
inside  the  surf  zone. 


I 


APPLICATION  OF  A NARROW  BAND  WAVE 
SPECTRUM  TO  THE  EXISTING  THEORIES 


The  expressions  just  considered  were  maae  as  general  as 
possible  in  order  to  afford  a simple  physical  understanding. 
Only  a monochromatic  wave  train  was  considered  end  only 
steady  state  solutions  were  allowed.  In  this  chapter,  un- 
steady terms  are  introduced  into  these  expressions  by  apply- 
ing linear  wave  theory  in  the  form  of  a simple  two  frequency 
component  wave  spectrum  simulating  a narrow  band  wave 
spectrum.  In  this  manner  an  unsteady  second  order  expression 
for  the  mean  water  level,  is  derived.  Next,  using  this 
simple  spectrum,  another  solution  for  q is  derived  using 
Longuet-H  i gg  i ns ’ application  of  the  Bernoulli  equation.  Then 
by  incorporating  this  spectrum  into  Iwagaki's  first  order 
solution  for  a sloping  bottom  LI972J  and  applying  the  Longuet 
Higgins'  approach,  a third  solution  for  q is  obtained. 
Finally,  the  application  of  the  spectrum  to  the  set  up 
phonnmonQn  inS!d6  Th0  SUPf  ZODS  I?  CnnSldSPGd* 

A.  WAVE  SET  DOWN 

Consider  a simplified  wave  spectrum,  specifically  consist 
ing  of  two  waves  which  have  nearly  tne  same  frequency  and 
wave  number  such  that  Oj  - a ^ ~ Aa,  and  kj  - k 2 - Ak  are  very 
small.  Also  for  simplicity  allow  tneir  respective  amplitudes 
to  be  identical  and  equal  to  — . This  ampli-^ude  is  chosen 

/7 

in  order  that  the  variance  of  the  simple  wave  spectrum  and 


the  variance  of  a monochromatic  wave  amplitude  are  the  same. 
The  wave  components  are  to  be  added  linearly  such  that  the 
particle  velocities  become  u = u.  + U2«  If  the  first  order 
velocity  potential  is  given  locally  by 

♦ = 22-  cosh  klhti).  i . I 2 

/?  ck  cosh  kh  ,1 


where  k.  is  used  to  indicate  direction,  i.e.  k /k  = cosa 

I X 

and  k /k  = sina,  then  in  general  u.  = - or 

y ^ I dx  . 


ag  I cosh  k(h+z) 


i k coshkh 

✓ 2 c 


cos ( k . X . - at ) . 


Since  u.  = u.  + u.  in  this  case,  u.  becomes 
' ' I ‘2  ' 


i,  cosn  k , ( h+z ) 

u.  = c— ' ! 

I /rr  k,  COShK.n 

/ C I I 


ccsCk.  X.  - G.t) 
1,1  I 


i^  cosh  k^Ch+z) 
~ coTrTTpi 


:os  ( X . X . - a_t ) . 
I . I 2 


Since  l^|  " 1^2  ~ subscripts  may  be 

dropped  and  (28)  can  be  written  as 


ag  r I cos  n k(  h-^z ) , , , 

u . = — — — cos(k.x.  - or 

I /2  c ^ cosih  Kh  If 

ki+Aki  GoshL  ( k + Ak)  ( h + z ) I]  , -n 

+ ( , 7Ti — ) rr-T — T ' T =^co s L ( k . +A  k . ) x . - ( 0+ Aa ) t j } . 

k +Ak  cosh(K-^AK;h  i i i 


li 


This  simplifies  to 


k A k 

2aq  i cosh  k(h  + z)  , i Aa  j. , ,, 

u.  = — 2—  ^ — cos(— X.-  t ) cos  ( k . X . -at ) . 

I ^ k cosn  kh  2 i z i i 


Similarly,  the  particle  velocity  in  the  z-direction,  w,  is 


obtained  and  is 


2ag  sinh  k(h+z)  , 

w = — 2 — cosC— X. 

— cosh  kh  z I 


/2  c 


Tj-  t ) s i n ( k . X . -at ) . 

2 I I 


These  expressions  stats  That  the  particle  velocities  are 
sinusoidal  with  a oeriodically  varying  amplitude.  The  modu- 
lating wave  is  of  low  frequency  and  with  a wavelength  that 
is  much  longer  than  those  of  the  original  two  base  waves. 

The  unsteady  velocities  specified  by  (35)  and  (36)  are 
substituted  into  (6)  and  the  radiation  stress  tensor  deter- 
mined. The  time  averaging  process  is  carried  out  over  the 
shorter  periods  of  the  base  waves.  Since  the  modulating 
wave  is  of  such  low  frequency,  its  effects  remain  in  the 
solution.  Specifically,  equation  (6)  becomes  term  bv  term 


c k . K . A k . X . , 

_ai:  ,11  Aa,,  , 

P U . ' U d Z = t — ^ C O S ( T t ) C O S ( 

I , = 2 2 


A k . X . 


t,  , 


I 2 

where  E = -y  pga  is  the  total  energy  densitv,  and  c and  c 
^ 9 

are  the  grouo  velocity  and  phase  velocity  respect i vs  I v ; 


c Ak. 

pdz  = pg(n  + h)’  + -I  - E(^  - 2)cos^(— ^ x.-  A? 


iT-  t ) ; 


v'33) 


and  the  remaining  unspecified  term  becomes 


M.M  . 

I J 


p ( n+h ) p ( n+h ) 


2 

c 


Substituting  (37),  (38),  and  (39)  into  (6),  the  radiation 

stress  tensor  becomes: 


S . . 

I J 


E 

c 


k.  k . 

I .] 
, 2 


Ak. 
cos  (^5-^ 


X . - 

I 


Aa 


t ) cos ( 


E(^  - 2) 

c 


- A k . 
cos  (— ^ 


Aa  ,, 
X . - t ) 

I 2 


( 40 ) 


, ^2  k. k . Ak. 

— !_  ccs(-jl. 

p ( r)  + h ) 0 k 


X . - 

I 


Aa 


Ak  . 


t ) cos  (• 


X : ” 


Aa  + . 


Since  only  second  order  effects  are  being  considered  and 
= 0(a'^)  the  last  term  may  be  neglected.  I is  obvious 
that  the  oroper  choice  of  assumptions  will  simplify  (40) 
immensely,  e.g.  consider  only  the  x-direction  components 
in  shallow  water. 


and 
a = 


XX 


k«  .i^Ak.  . n 

= LE(^)^+  £j  cos^(^  ^i"  — ^ f ' 


f the  waves  aporoach  the  shoreline  such  that  The  angl 
),  this  expression  further  reduces  to 


S 

XX 


c 

+ cos(Ak.x.  - Aat)  + . 

I i 2 


(4  , ) 


This  last  expression  states  that  the  radiation  stress,  in 
this  circumstance,  is  composed  of  a steady  stress  term  and 
of  one  That  is  periodic. 

This  derived  form  of  the  radiation  stress  is  now  used 
to  find  a solution  for  n.  Recalling  equations  (10)  and 
(II),  and  as  before  neglecting  the  friction  term  and  assum- 
ing that  the  current  gradients  in  the  x-direcrion  are  small, 
the  x-component  of  these  expressions  oecome: 


3x 


P IT 


3y  cross  differenriating.  The  M term  can  be  eliminated  and 

X 

" p rn  ^ 17^"  IT-  - - 7-r-  ^^2) 

<3  1 ax 


is  obtained.  Assuming  a flat  bottom  so  that  h is  considered 


constant,  (42)  becomes 


2-  T_ 

3 n ^ ^ 3^ri 

- p — ^ + ogh  — ^ = 


This  is  obviously  the  long  wave  equation  which  is  being 


*orced  by  the  radiation  stress  Term. 


nas  a soluTion  in 


I 


-ne  homogeneous  case  of 


r)_  = a cos(Ak  x - Apt)  . 


For  a particular  solution,  make  the  assumption  that  led  to 


4- 


equation  (41),  i.e.,  consider  shallow  water  and  a = 0.  The 

assume  a particular  solution  of  the  form 

r)  = A cosCAk  x - Aot)+  B sinCAk  x - Act)  . 

p X X 

By  substituting  -l-his  into  equation  (43),  the  solution  for  A 
and  B are  obta i ned 

EAk  ^ 

A = ^ , 

p(Aa~  - ghAk^) 

3=0. 


The  solution  for  n is  then 


Ca  - 


2^9 


(^  - gh) 
Ak^ 


-H  cos(Ak  X - Aat)  , 

X ' 


( 4.i  ) 


which  is  a wave  with  an  amplitude  that  is  proporrional  to 
the  steady  state  solution  of  Long uet-H  i gg  i ns  and  Stewart, 
equation  (15). 

If  the  same  assumptions  as  before  are  made  except  tha- 
the  bottom  Is  allowed  to  vary  as  h = mx  wnere  m is  some 
constant  slope,  a more  difficult  croolem  is  encountered. 
Equation  (42)  becomes 


3n 

•P  4 + Pgm  -5- 

^t^  3x 


pgmx 


0 _ 

S'n  _ 
3x" 


XX 


(45) 


3x‘ 


This  equation  again  has  the  characterisrics  cf  the  long  wav 
equation  that  is  being  forced  by  the  radiation  stress.  An 


analytical  solution  of  this  equation  would  provide  a general 
form  for  fi ; however,  since  the  depth  is  al  lowed  to  vary  with 
X,  the  expression  for  the  radiation  stress  term  becomes  com- 
plicated and  obtaining  a particular  solution  for  (45)  becomes 
difficult. 


B.  WAVE  SET  DOWN  USING  THE  BERNOULLI  INTEGRAL 

Recalling  that  the  second  method  of  Longuet-H i gg i ns  to 
determine  a solution  for  n resulted  in  the  ecuation 


A-  I r.-  2 ^ - 2 -2,  -,1 

- 23  t^^^x  ^ ^y  - '''  ^z  = 0"2  ' 


a second  unsteady  solution  may  be  obtained  by  applying  the 
simple  wave  spectrum.  From  equations  (35)  and  (26)  u 

X 

u ^ , and  w become 


o 22  1,2,,,..  _Ak  . 

-2  a g , X.2  cosh  k(h■^2)  2,  x Aa  ,, 

2 k k 2 2’ 

c c o s h k h 


n 2 = (^)2  cosh2k(h  + 2)  ^032^:^  y - t)  , 

\/  2 k ^ 9 ' / ' 


cosh  kh 


n 22  -1.2,  , u,  N oAk  Ak  , 

-2  ag  s I nh  k ( h+z ) 2,  x . v Aa., 

w = — H 5 cos  (-^  X +-y-  y - _ t)  , 

c cosh^kh 


I -2-2-2 
and  - ■^(  u +u  -w  ) ^ 

2g  X y z=0 


22  x_  -jAk 
3gr,x,2  2,  X Aa^, 

— ^{(-j^)  cos  (—rr~  x~  ~ t) 

2c"  " " 

k _ A k , 

+ (-^)2c03“(-t^  V-  ^ t) 


- tanh  kh  cos' 


!■  X ^ y Aa  , , , 

( — s—  X + V - -s-  t ) } 


Again  assuming  a = 0 so  that  = l»  snd  k^/k  = 0 and 

2 

recalling  that  c = g/k  tanh  kh,  (46)  becomes 


n = - 4 ^ ^ ' Cl  + cosCAk  X - Aat)]  . (47) 

2sinh2kh  x 


This  solution  is  identical  to  the  solution  of  Longuet-Hig- 
gins  (25)  except  for  the  unsteady  term.  The  set  down 
fluctuates  at  Aa  = Oj  - 02',  which  is  the  difference  of  the 
radial  frequencies  of  the  two  component  wave  spectrum.  This 
low  frequency  oscillation  is  often  described  as  surf  beat. 

The  amplitude  of  n varies  from  twice  the  amplitude  of 
(25)  when  the  two  base  waves  are  in  phase,  to  zero  when  the 
base  waves  are  out  of  phase.  Hence  the  set  down  can  be  con 
sidered  the  superposition  of  a steady  component  and  an 
unsteady  component. 


C.  WAVE  SET  DOWN  ON  A SLOPING  BEACH 

In  an  effort  to  more  accurately  represent  shoaling  wave 
transformations,  Iwagaki  lI'572C  implicitly  considers  con- 
stant bottom  slope  resulting  in  solutions  to  the  long  wave 
equation  involving  the  Bessel  functions.  He  obtained  first 
order  solutions  for  n and  u in  terms  of  the  asymptotic  form 
of  the  Bessel  and  Neumann  functions.  It  seems  reasonable 
to  assume  that  by  incorporating  these  solutions  into  the 
Bernoulli  integral,  a second  order  solution  for  fi  can  be  ob 
tained  which  includes  a mere  realistic  representation  of  th 
bottom  effects.  Assuming  shallow  water  conditions  with  the 
deoth  given  by  h = mx  and  progressive  waves. 


■ ^1/^ 


u 


[j.(X)sin  Jt  + N.(x)cos  otj 


; 48  ) 


n 


(49) 


= — CJ  (x)cos  at  + N (x)sin  at'J  , 

/2  ° ° 

2 

2 ax 

where  x ~ ^ • The  asymptotic  expansions  of  the  Besse 

and  Neumann  functions  are: 


f TTO) 


J (Oi  ) ^ •*/■=—  COS  (w 
n 


N (o)) 
n 


~ V-  sin 

f Trcj 


( 0)  - 


nTT 

2 


niT 

2 


T> 


T> 


( 50 


Substituting  these  expressions  into  (48)  and  simplifying 
yields 


= 

/2 


Cs  i n (X  - -r 


3Tr 

4 


at ) J 


(51) 


and  from  the  kinematic  free  surface  boundary  condition 


_ 9n 


at)j 


(52) 


Y *1 

Introducing  the  wave  spectrum  as 


averagi ng 

the 

res  u 

It, 

produces 

II 

o 

II 

N 

CM 

X 

I Z3 

2 

£_£ 

2mx 

/gmx 

2 

TTO 

L( 

,2Aax 

r COS  ( 

/gTnx 

and 

-2 

w 

2 2 , 

a a /qmx 

2xTTa 

[1 

,2Aax 
+ cos( 

before. 


+ Aat)Il 


+ Aat)3 


squaring  and  time 


(53) 


( 54) 


By  assuming  that  ri2  's  in  infinitely  deep  water  as  before, 
tne  second  order  solution  for  ii  is 


[_2_  - a^nCi+cosC^  + Act)]  . (55) 
2g7rax  mx  gmx 


D.  WAVE  SET  UP 

Because  the  nature  of  the  solution  for  n using  linear 
theory  was  the  superposition  of  a steady  state  component  and 
an  unsteady  component,  it  is  reasonable  to  assume  that  this 
condition  persists  across  tne  breaker  line  and  the  set  up 
resulting  in  the  surf  zone  from  the  wave  spectrum  has  a 
similar  makeup.  The  steady  state  component  is  then  derived 
separately  in  the  same  manner  as  in  Chapter  II  with  the 
radiation  stress  being  now  defined  as  the  steady  state  por- 
tion of  equation  (41),  or 


XX 


( 56) 


and  aga  i n 


n.  = - kh  + C , 


( 57) 


where  this  time  the  proportionality  constant  is  again 


k;  = 


3 Y 


Since  Hu  = (57)  becomes 

J 10  0 


= (h,-h)+nu 

s b b 


58) 


d 


The  unsteady  component  is  assumed  to  be  periodic  in 
character;  however,  unlike  the  higher  frequency  waves  of 
the  train,  observations  suggest  that  it  is  not  attenuated 


as  it  approaches  the  shoreline  but  is  reflected  to  some 
extent.  The  formulation  is  then  that  of  the  long  wave 
equation  as  given  by  Stoker  Cl966j, 


2-  2 

5 ^ - gh  ^ ^ 

8t^  9x^ 


gh 


9 h 9ri 
9x  9x 


0 , 


( 59) 


and  the  solution  is  given  by  Guza  and  Bowen  Cl977j  as 


n = a{J  (x)sin  Aat  + Nl  (y)cos  Aot  + 
o o ' 


f 

■» 

i 

i 

i 


'‘CJQ(X^sin(Aat-e)-N^(X)cos(Act-£)J}  (60) 


where  a is  the  amplitude  of  the  incoming  wave  (in  this  case 
3 = ^ reflection  coefficient,  and  £ a phase 

shift  in  the  reflected  wave.  The  frequency  of  this  wave 
must  be  the  same  as  that  of  the  modulating  wave  derived  cur- 
side  the  surf  zone  or  Ac.  By  substituting  +he  asymptotic 
form  of  the  Bessel  and  Neumann  functions, 

{cos ( x-Aat  - j) + r cos ( x-Aat-0 ) } , (61) 

where  9 = -^  +c.  Since  the  set  up  shoreward  was  assumed  to  be 

4 

a composite  of  the  two  solu'l'ions  (58)  and  (61), 


= k ( h ^-h  ) -rn , 


{I 


•^Lccs  (x-Acrt-  -j) -i-rcos  ( x-Aat-6  ) 2 } 


(62) 


IV.  COMPARISONS 


The  results  obtained  here  are  compared  with  the  labora- 
tory experiments  performed  by  Bowen,  et  al.Cl968];  Figure  2. 
They  obtained  experimental  results  from  a controlled  wave 
tank  experiment,  and  in  each  case  only  a monochromatic  wave 
was  considered.  Their  results  compare  quite  favorably  to 
the  steady  state  solutions  of  Longuet-H i gg i ns  and  Stewart 
Cl  962]. 

Since  it  was  assumed  that  The  variance  of  the  simple 
wave  spectrum  was  equal  to  the  variance  of  the  monochromatic 
wave,  the  derived  expressions  differ  from  existing  theory  by 
that  amount  which  is  contributed  by  the  moaulating  wave  of 
the  unsteady-  portion.  Comparing  Longuet-H i gg i ns ' solution 


2 s i n h 2 kh 


with  the  unsteady  solution  (47) 


I r,  , ,, 

" - - I sinh  2kn-'  * - iatlj  , 


it  is  found  that,  by  averaging  tne  effect  of  the  unsTeady 
component,  they  are  the  identical.  As  mentioned  earlier, 
the  total  effect  of  the  unsteady  motion  is,  in  this  case, 
to  produce  a set  down  that  oscillates  periodically  from 
twice  the  steady  solution  to  zero  times  it,  i.e.,  no  set 


down  at  all.  Thus  a group  of  waves  which  are  close  in  fre- 
quency and  wavelength  may  be  said  to  produce  in  the  near 
shore  region  a fluctuation  in  the  mean  sea  level  seaward  of 
the  breaker  line.  This  is  also  demonstrated  by  the  unsteady 
solution  obtained  by  allowing  for  a sloping  bottom, eq uat i on 
(55).  It  too  consists  of  a steady  component  and  an  oscillat- 
i ng  component. 

The  initial  solution,  equation  (44),  on  the  other  hand 
is  comprised  of  only  an  oscillating  component.  However,  this 
solution  was  derived  in  a simplistic  manner  by  assuming  That 
the  depth  remained  constant  and  by  solving  the  resulting 
differential  equation.  While  this  demonstrates  that  the 
mean  sea  level  oscillates  in  a wave-like  fashion  when  driven 
by  a group  of  waves,  it  does  little  to  provide  insight  into 
the  near  shore  region. 

The  set  down  can  be  expressed  in  shallow  water  in  terms 
of  deep  water  conditions  by  using  conservation  of  energy 
flux. 


E c = pga  c = constant  , 
9 - 9 


w h e e 


pga 


o 2o 


rr  pga  . /g  n 


° i I 


or 


/ 


The  subscript  o refers  to  deep  water  and  the  subscript  s 
refers  to  shal low  water.  The  set  down  in  shal low  water 
by  equation  (25)  is  then  given  by 


2 1/2,  ,-3/2 

a g ( mx ) 


where  in  shallow  water  sinh  2kh  -►  2kh  and  h = mx  for  a con- 
stant sloping  bottom. 

The  set  down  derived  from  linear  theory  can  be  compared 
to  that  derived  from  linear  theory  for  a sloping  bottom, 
equation  (55),  by  also  expressing  the  shallow  water  amplitude 
for  sloping  bottom  solutions  in  terms  of  deep  water  condi- 
tions. This  is  accomplished  by  patching  Stokes  solution 
at  the  off  shore  to  the  shal low  ater  solution.  For  a smooth 
match  it  is  required  that,  from  Friedrich  Cl948ll, 


2 

a •rr 
o 

2m  • 


The  steady  state  solution  of  (55)  can  now  be  expressed  In 
terms  of  deep  water  conditions 


2 2 
- g'/^(mx)"^^^  + -2 — (gmx)' 

f)  I + —2 — ( gmx  ) “ ' '^^  . 


Hence,  the  sloping  bottom  solution  results  in  slightly 
greater  set  down  as  compared  with  linear  *heorv.  Numerics 


results  of  these  equations,  (63)  and  (64),  are  plotted  in 
Figure  2 . 

Close  to  the  break  point,  most  theoretical  results  fail 
to  compare  favorably  with  the  experimental  results.  Theory 
increases  set  down  rapidly  near  the  break  point  while  ex- 
perimental results  tend  to  flatten  out.  This  is  due  most 
likely  to  the  failure  of  linear  theory  to  adequately  describe 
breaking  waves. 

Although  simplistic  in  form,  the  application  of  the  two 
frequency  wave  spectrum  demonstrated  that  fluctuating 
values  are  obtained  for  the  set  down  and  set  up  phenomena, 
i.e.,  time  dependent  solutions  were  obtained  for  these  con- 
cepts. Hence,  a group  of  waves,  similar  In  frequency  can 
be  expected  to  produce  a periodic  variation  in  the  mean  sea 
level  wnich  becomes  most  apparent  when  they  encounter  a 
beach  . 


38 
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V.  CONCLUS I ONS 

In  this  study,  unsteady  solutions  for  set  down  and  set 
up  were  derived  using  a simple  two  component  wave  spectrum. 
The  spectral  wave  components  were  described  using  linear 
wave  theory  for  both  horizon+al  and  sloping  bottoms.  The 
sloping  bottom  solution  involving  the  Bessel  functions  gave 
a slightly  greater  set  down  comoareG  with  mat  of  rne 
linear  wave  theory  solution.  However,  soth  unsteaay  solu- 
tions for  set  down  showed,  that  to  at  least  a first  aoproxi- 
mation,  a steady  and  a fluctuating  comoonent,  the  steady 
component  being  identical  to  the  earlier  steady  state 
sol u 1 1 on  s . 

The  set  down  at  the  breaker  I i ne  acts  as  the  boundary 
condition  driving  the  se+  up  inside  the  surf  zone.  The  set 
down  solution  showed  that  the  steady  and  unsteady  components 
are  simply  additive.  Hence,  inside  the  surf  zone  a solu- 
tion is  composed  of  the  steady  set  up  ana  a long  wave,  i.e., 
sjrf  beat,  which  is  driven  by  tne  fluctuating  condition  at 
the  brea  ker  line. 
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